Trace formula for dielectric cavities : I. General properties 
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The construction of the trace formula for open dielectric cavities is examined in detail. Using 
the Krein formula it is shown that the sum over cavity resonances can be written as a sum over 
classical periodic orbits for the motion inside the cavity. The contribution of each periodic orbit 
is the product of the two factors. The first is the same as in the standard trace formula and the 
second is connected with the product of reflection coefficients for all points of reflection with the 
cavity boundary. Two asymptotic terms of the smooth resonance counting function related with 
the area and the perimeter of a convex cavity are derived. The coefficient of the perimeter term 
differs from the one for closed cavities due to unusual high-energy asymptotics of the S-matrix for 
the scattering on the cavity. Corrections to the leading semi-classical formula are briefly discussed. 
Obtained formulas agree well with numerical calculations for circular dielectric cavities. 

PACS numbers: 42.55.Sa, 05.45.Mt, 03.65.Sq 



I. INTRODUCTION 

Only a limited number of quantum models permits ex- 
act analytical solutions. All others require either nu- 
merical or approximate approaches. One of the most 
useful approximation and widely used for non-integrable 
multi-dimensional quantum systems is the semiclassical 
one based on different types of trace formulas developed 
in the early 70's [ll,i,[lQ. 

A general trace formula relates two different objects. 
Its left-hand side is the density of discrete energy levels, 
Ejn, of a quantum system 
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The right-hand side of the trace formula is the sum of 
two contributions 



i— classical 



{E) ^ d{E) + d^°"'\E) 



(2) 



The first term, d{E), is the smooth part of level density 
given by a series of the Weyl terms (see e.g. pj] and 
references therein). By definition, d{E) = dN{E)/dE 
where N(E) is the mean number of levels with energies 
En < E called also the smooth counting function. For 
two-dimensional billiards leading contributions to N{E) 
are the following 
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where k = \/E is the momentum, A is the area of the 
billiard cavity, and L is its perimeter, r in this formula is 
related with chosen boundary conditions. For example, 
r = ±1 for, respectively, the Neumann and the Dirichlet 
conditions. 



The second term in ^, d^""^^ (_B) , is the oscillating 
(or the fluctuating) part of the level density given in the 
leading order by the sum over classical periodic orbits 
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periodic or 



bits 



where dp{E) is the contribution of a given periodic orbit 



dp{E) = 



„iSp(B)/ri 



(5) 



Here Sp{E) is the classical action along the periodic or- 
bit (for billiards Sp{E) = kip with Ip being the length 
of the periodic orbit) and the prefactor Cp can be com- 
puted from classical characteristics of the given periodic 
orbit (see below (f24| and ([25| ). The trace formula is the 
statement that 
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In an overwhelming majority of cases the trace formulas 
were applied to closed quantum systems with discrete 
energy levels. Though it is known (see e.g. Q) that the 
trace formula and its modifications can also be applied for 
open quantum systems, only a small number of examples 
has been considered so far 0]. 

This is the first of a series of papers whose purpose is 
to demonstrate the usefulness of the application of trace 
formulas to open quantum systems, in particular, to di- 
electric cavities where a bounded domain is filled with 
a dielectric material with refractive index n > 1 while 
the exterior is a media with refractive index 1. The in- 
creasing interest to the investigation of such models to 
a large extend is related with their potential importance 
as resonators in micro-lasers (7|, |8| and the experimen- 
tal possibilities of observing geometrical characteristics 
of the cavity from their lasing spectra [3, [lO] ■ The con- 
struction of trace formulas for dielectric cavities has been 
briefly discussed in 10]. We also mention that for open 
chaotic resonators a variant of a trace formula has been 
developed in [ll| . 



Open systems, in general, have no discrete spectrum. 
So it is not immediately clear what should be the left- 
hand side (the spectral part) of the trace formula. In Sec- 
tion|TT]we discuss the Krein formula [12] which is the main 
theoretical tool for open systems. This formula relates 
the excess density of states for an open cavity with the 
derivative over energy of the determinant of the S-matrix 
for the scattering on the cavity. From the Krein formula 
it follows that the spectral part of the trace formula con- 
sists of the sum over all complex poles of the S-matrix 
widely called resonances or quasi-stationary states. 

The right-hand side of the trace formula should include 
a smooth part and a sum over periodic orbits (cf. ^). 
The general formalism like the multiple scattering ap- 
proach, used by Balian and Bloch in their construction 
of the trace formula for closed cavities [ij , is not yet fully 
developed for dielectric cavities. The main difference of 
dielectric cavities from closed ones is that for the later 
one knows boundary conditions on the cavity boundary 
but for the former boundary conditions are fixed only at 
infinity. Due to this fact the resulting integral equation 
(I13p includes the integration over the cavity volume and 
not over its boundary as for the closed case which compli- 
cates the general construction. In Section Hill the form of 
leading periodic orbit contributions is fixed from physical 
considerations. 

To get more insight to this problem, in Section IIVI 
a simple example of integrable circular dielectric cav- 
ity is considered in detail. Despite the calculations are 
straightforward, they permit us to calculate in Section FVl 
the smooth part of the counting resonance function and 
in Section IVII the periodic orbit contributions. In Sec- 
tion |Vlll we consider main corrections to the asymptotic 
results and, in particular, demonstrate how the Goos- 
Hanchen shift [lj| manifests in the trace formula. All 
obtained formulas agree well with numerical calculations. 

Though these results were obtained, strictly speaking, 
only for circular dielectric cavities, all main steps leading 
to them are quite general and we conjecture that they re- 
main valid for dielectric cavities of arbitrary shape. The 
detailed comparison of the derived trace formulas with 
numerical calculations for cavities of different shapes is 
postponed for a future pubhcation [14 |. 



II. KREIN FORMULA 

Throughout this paper we focus on a planar two- 
dimensional domain V. The wave equations widely used 
for such dielectric cavities are the following (see e.g. |15| ) 






when x £ V , 
when x ^ V . 



(7) 



In electromagnetism we rather consider the wavenum- 
ber k = vE, which is related to the frequency uj of the 
wave through k = co/c, c being the speed of light in the 
vacuum. For simplicity we shall suppose that the wave 



function, ^(x), and its normal derivative are continu- 
ous along the boundary of V. In electrodynamics these 
boundary conditions correspond to transverse magnetic 
(TM) modes inside an infinite dielectric cylinder with 
cross section V [l5J . Transverse electric (TE) modes may 
be treated similarly and will be considered elsewhere. 

We stress that Eqs. ([7]) describe correctly electromag- 
netic fields only for an infinite dielectric cylinder. For 
a real 3-dimensional cavity with a small height, /i, they 
have to be modified. The authors are not aware of the 
full 3-dimensional treatment of such problems. The usual 
approach consists to consider the refractive index, n, in 
these equations not as a constant but as an effective re- 
fraction index, n = neff{kh), for the motion inside a 
2 dimensional dielectric slab (see e.g. 16] and refer- 
ences therein) but errors of a such approximation are 
not known at present. 

The possibility of using such cavities as resonators is 
related with the phenomenon of the total internal reflec- 
tion. It is known (see e.g. [l5|) that if the TM plane wave 
inside the cavity is reflected from a straight boundary, the 
(Fresnel) coefficient of reflection has the following form 
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when \9\ < 9r 



, when \9\ > 9c 



(8) 

Here 9 is the angle between the direction of the incoming 
wave and the normal to the boundary and the critical 
angle 



9c = arcsin 1/n 



(9) 



When 16*1 > 9c |^(^)| — 1 and the wave is completely 
reflected from the boundary which may lead to formation 
of long-lived states. We mention that this expression 
becomes less efficient for curved boundary close to the 
critical angle (see e.g. pL7|). 

Eq. ([7]) admits only continuous spectrum and its prop- 
erly renormalized eigenvalue density, d{E), has to be a 
smooth function of energy, contrary to closed systems 
where the level density is a sum of delta's peaks (cf. H])). 
It is convenient to rewrite ([7]) as the Schrodinger equation 
with a potential 

{A + E)-i'{x)^-XV{x)'^{x) (10) 

where the coupling constant 

X^E{n^-l) (11) 

and the potential V{x) is non-zero only inside the cavity 



V{x) 
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(12) 



Except of unusual dependence of the coupling constant 
on the energy this equation describes the motion in a 



finite-range potential so standard methods (see e.g. |18| ) 
can be applied to analyze it. In particular, its solution 
corresponding to the continuous spectrum is defined by 
the following integral equation (see also [19J ) 



^(x) 
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(13) 



where k is the momentum vector of the incoming wave 
with coordinates fc(cos 6',sin 6), Go(A:|a; — y\) is the free 
Green function of the left-hand side of ([T(I| 



Go(z) = ii/^')(z), 



(14) 



r(i) 



and Hq (z) is the Hankel function of the first order. 

This equation may serve as a starting point of the mul- 
tiple scattering method similar to the one discussed in [l| . 
The important difference between these two cases is that 
the integration in (|13p is performed over the whole cavity 
volume which complicates the iteration procedure. Nev- 
ertheless, Eq. (fT5|) has the standard form of the Fredholm 
equation and its solutions are well defined similar to the 
usual case when coupling constant A is independent of en- 
ergy. In particular, from (jlSp it follows that the S-matrix 
for the scattering on the cavity has the form 

s{e,e') = d{0-e') + ^f{9,9') (15) 

47r 

where 6 and 9' are the angles determining the directions 
of, respectively, incoming and outcoming waves and 



f{d,d') 



-'^■'%(y )dy . 



(16) 



Here k' is the momentum of the outgoing wave with co- 
ordinates k{cos9' ,sm9'). 

The importance of the S-matrix lies in the fact that the 
excess density of states for open systems can conveniently 
be written through it by using the Krein formula (see |12| 
and [2y, |2l|). This formula relates the density of states 
of two operators: the first with a short-range potential, 
d{E), and the second without it, do{E). In physical terms 
it reads 

d{E) - do{E) = ^^ In det S(i?) (17) 

where det S(E) is the determinant of the S-matrix for the 
scattering on the potential. For clarity in Appendix [X] a 
physical 'derivation' of this formula is given in the sim- 
plest case of one-dimensional systems. 

It is also known (see e.g. [l^, sect. 12, and [23. |23|') 
that function det S{E) for the scattering on a finite-range 
potential is a meromorphic function in the complex k = 
yE plane with (in even-dimensional spaces) a cut along 
the negative /c-axis 



det S(£;) 
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(18) 



where km with Im fcm < denote the positions of the 
poles of the S-matrix on the second sheet of energy sur- 
face and the product is taken over all such poles. (Terms 
which ensure the convergence of this infinite product are 
not written explicitly). Here /(fc) is a certain function 
without singularities in the cut complex plane related 
with the asymptotics of the S-matrix when E —^ oo. Due 
to the symmetry k — > — fc of ([7]), if km is a pole, — fcj^ will 
also be a pole which is implicitly assumed in ^TE\i . 

Therefore, the left-hand side of p?)) can be written as 
the sum over all poles of the S-matrix 

m 

where E„i = k^ and g{E) = /'(fc)/27r. In mathematical 
literature such formulas are known as Poisson formulas 
(see [23 and references therein). 

From (fTg]) it is clear that poles of the S-matrix, espe- 
cially those close to the real axis, play an important role 
in d{E) as they produce peaks in the otherwise smooth 
background. 

The positions of these poles, called also resonances or 
quasi-stationary states, are in principle determined from 
Eq. ([7|) by an analytic continuation of the S-matrix from 
the real axis to the complex plane, in a similar way as 
in [24] . Another way is to impose the outgoing radiation 
conditions (see e.g. pJi]) 
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(20) 



These conditions select a well defined set of complex 
eigenvalues Em whose imaginary parts correspond to the 



resonance lifetime Tm, Im E„ 



-l/2r„ 



These arguments demonstrate that the spectral part 
of the trace formula (the analog of H])) may be written 
as the sum over all poles (= resonances) of the S-matrix 
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(E - Re E„ 
Im km 



(Im Em)^ 
(21) 



When Im i?„ 
tribution. 



(fc - Re fc„)2 + (Im km)^ ■ 

one recovers the usual (5-function con- 



III. GENERAL PROPERTIES OF TRACE 
FORMULAS FOR DIELECTRIC CAVITIES 

To find the right-hand side of trace formula one has 
to express it through the Green function. In this section 
we will only consider the real continuous spectrum of 
Eq. ([7]) and eigenfunctions associated with it. We will 
not deal with functions associated with the resonances 
defined by boundary condition (|20p . A minor difference 



with the usual case of closed systems consists in the fact 
that the eigenfunctions are orthogonal not with respect 
to the standard scalar product / ^^^ {r)'^E2 (r )dr but to 
a modified one 



\f)^*E,{r)^E"{rW= 6{E' - E") 



(22) 



where the integration is performed over the whole space. 
Here n{r) = n for points inside the cavity and n(f) = 1 
outside it. This relation can easily be checked from the 
main equation Q. 

This modification leads to the following formal relation 
between the density of states and the trace of the Green 
function (cf. e.g. [ll]) 



d{E) = -- I n^{r)lTii GE{r,r)dr 



(23) 



where the integration is extended to the whole space. 

Exactly as for closed systems (see e.g. [ll, [2, S|)j the 
dominant contribution to the fluctuating part of the den- 
sity of states comes from classical periodic orbits. When 
saddle point is considered, it is clear that the only in- 
formation that a trajectory may have about a boundary 
is contained in the reflection coefficient from this bound- 
ary. Therefore, the principal difference with the closed 
systems is that the contribution of a given periodic orbit 
has to be multiplied by the product of reflection coeffi- 
cients ([S|) for all bounces with the cavity boundary (which 
we called the total reflection coefficient) and the later can 
be less than unity. 

These simple arguments demonstrate that for dielectric 
cavities the dominant contribution of a periodic orbit to 
the density of states when k -^ oo has the following form 

• For an isolated primitive periodic orbit p repeated 
r times 

where Ip, Mp, ^p, Rp are, respectively, the length, 
the monodromy matrix, the Maslov index, and the 
total reflexion coefficient for the chosen primitive 
periodic orbit. 

• For a primitive periodic orbit family 



dp{E) 
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where Ap is the area covered by periodic orbit fam- 
ily, (Rp) is the mean over family value of the total 
reflection coefficient, and all other notations are as 
above. 

The dependence on n of the prefactors in these for- 
mulas is related with the fact that inside the cavity the 
momentum equals nk. 



Now we have all ingredients of the trace formula for 
dielectric cavities except the smooth (Weyl) terms ([3|). 
It is clear that its value depends on how many S-matrix 
resonances are included in the right-hand side of the trace 
formula ((2T|) . Energy eigenvalues of closed systems are 
real and they all have to be included. For open systems 
resonance energies are complex and a natural approach 
is to take into account only resonances whose imaginary 
part is restricted, e.g. — Im km < 7. Other resonances 
(if any) have to be considered as a smooth background. 
This separation of poles is, to a large extent, arbitrary 
which manifests in the fact that the smooth part of the 
resonance density will be now a function of 7 whose cal- 
culation is a difficult problem. 

For 2-dimensional open chaotic billiard problems with 
holes there exist strong arguments [25,, 2§l that the lead- 
ing term of the Weyl law is non-trivial 



N{E) - C(7)fc'' 



(26) 



where ;y < 2 is related with the fractal dimension of the 
trapped set of classical orbits. 

For open dielectric cavities arguments lea ding to ([26]) 
can not be directly applied and it appears [27| that in 
this case standard estimates ([3]) with different constants 
are valid. In general, there exist two types of resonances 
whose behavior is different when the system is somehow 
" closed" . The first (sometimes called Feshbach or inter- 
nal resonances) are resonances which tends to eigenval- 
ues of the corresponding closed system and the second 
ones (called shape or outer resonances) are those whose 
imaginary part remains non-zero when a system is closed. 
The both types of resonances exist for dielectric cavities 



(see e.g. [16|). In all cases with convex shaped cavities 



considered by us these two groups of resonances are well 
separated and there exists a clearly defined value of 7max 
such that all first type resonances have — Im km < 7inax- 
It does not contradict [2J] where a Fractal Weyl law was 
argued to be a good description of a part of the resonance 
spectrum with small imaginary part. In the next section 
we shall argue that when all these resonances are taken 
into account the smooth counting function N{E) defined 
as the mean number of resonances with Re En < E = k"^ 
and — Im fc,„ < 7,nax is similar to ([3]) but with following 
modifications 

N{E)^^+f{n)^ + 0{l) (27) 

where A and L are, as above, the area and the perimeter 
of the cavity but f{n) is 



f{n) = 1 + 
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^R{t) - - 



dt 



f^ + l 



and 
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Rit) 
(28) 



(29) 



The function r{n) can be expressed through eUiptic in- 
tegrals. It is monotonic function of n starting from 1 for 
n — 1 and tending to n(4/7r — 1) for large n. 

Finally, the trace formula for dielectric cavities has the 
following form 
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n 



E 



Im kri 
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2k{d{E) + Y,KiE)+d;{E)]} 



(30) 



Here dp{E) for isolated periodic orbits is given by 
and for periodic orbits from a family it is defined in (j25p . 
The factor 2k in the right-hand side is introduced as we 
found it more convenient to work with the density of 
states as a function of momentum. 

To see the effect of periodic orbit it is usual to mul- 
tiply the both sides of the trace formula ([30]) by a test 
function, e.g. by e"'*^"', and integrate over fc in a certain 
window ki < k < k2 which includes many resonances. 
The dominant contribution in the left-hand side of this 
formula is the sum over resonances whose momentum is 
restricted fci < Re fc < fc2 and < — Im km < 7max 



E< 
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(31) 



where the summation is done over periodic orbits, p, and 
it is implicitly assumed that term with p = corresponds 
to the smooth part of the trace formula. Here Ip de- 
notes the integral the individual periodic orbit contribu- 
tion given by ([24]) or ^ 
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dp{E)2kdk 



(32) 



fci 



which is strongly peaked at the periodic orbits length. 
Eq. ([3T|) is called the length density. 

Periodic orbits are not the only contributions to the 
trace formula. For polygonal cavities important contri- 
butions are given by diffractive orbits which go through 
singularities of the boundary (see e.g. [28J and references 
therein) . Usually their individual contribution is smaller 
than those of periodic orbits. The careful calculation of 
such corrections requires the knowledge of the diffraction 
coefficient on dielectric singularities which is not avail- 
able analytically. We shall also not discuss here creeping 
orbits (see e.g. [2^ and references therein) corresponding 
to the external motion along cavity boundary as they are 
responsible only for shape resonances and their contribu- 
tions are small. 



IV. CIRCULAR CAVITY 

The circular dielectric cavity is one example of a two- 
dimensional cavity with an explicit analytical solution 
and it is instructive to illustrate the above general for- 
mulas in this simple case. 



The Green function G{'r,f',E) for this problem has 
been written e.g. in [la |. From that formulas it follows 
that inside the circular cavity with radius R the Green 
function G{f, r, E) with r < R has the form 
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where 

9i7i{r) = Jm{nx)H^\nkr) - J,n{nkr)H^^ (nx) 
and outside it (when r > R) it is 
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with 



9out{r) = Hl^Xxyjmikr) - Jm{x)H';^\kr) 



(36) 



(1)^ 



Here x = kR, Jm{z) (resp. Hm [z)) stands for the Bessel 
function (resp. the Hankel function of the first kind), and 



SmK^) = n-—[nx) 

Jrn H, 
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(37) 



Quantities Smix) are of special importance as the posi- 
tions of resonances, Xm = kmR, are determined by their 
complex zeros (see e.g. |16|) 



^m\'^m } U 
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Using the formula |3C 



A^{t)B^{t)tdt 
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-A,{x)B,{x) (39) 



valid for any two Bessel functions A^{z) and Bi,{z) of the 
same arguments and performing straightforward calcula- 
tions one gets 
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where d{E) is formally 



d{E) = 27r 



Im G{f,r)rdr (41) 
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FIG. 1: a) Resonance spectrum for the dielectric circular cav- 
ity with n = 1.5 and Re kR < 100. b) The same but for n = 2 
and Re fci? < 40. 



and 



do{E) = 2tt I Im Go(r,r)rdr 



with the free Green function 



G^{r,r') = ^Hl^\k\r-r'\ 



(42) 



(43) 



In the above formulas it is assumed that the cut-off A ^ 
R and the limit A ^ oo is taken after the subtraction 
d{E)-do{E). 

The calculation of the S-matrix for the circular dielec- 
tric cavity is also well known [3l| and one finds 



^m[E) 7TT— — 



Hm {x)Sm{x) 



(44) 



with s,ra(x) are defined in (|37p . It is straightforward to 
check that PO)) can be written in the form 



d(E) - do{E) 



J_d_ 
2m dE 



In dot S(£;) 



(45) 



which is the Krein formula (I17p for the dielectric disk. 

In Figs. [1^) and b) we present numerically computed 
from (|38p positions of complex momenta of quasi-bound 
states for dielectric circular cavities with n = 1.5 and 
n = 2. Clearly, resonances in this case are organized in 
families corresponding to conserved quantum numbers of 
radial and angular momenta (see e.g. [16]). Let notice 
that for large Re kR imaginary parts of all resonances 
obey the inequality 



Im kR < 7,] 



where 



1 



7max(n-) = —In - 

2n n — 1 



(46) 



(47) 



7max(l-5) ~ .5365 and 7niax(2) ~ .2765 in a very good 
agreement with Fig.[T^) and b). This bound corresponds 
to the imaginary part of zeros of Sm{x) with m = and 



Re('kR) 



FIG. 2: Outer (shape) resonances for the dielectric circular 
cavity with n = 1.5 Dashed line indicates the lowest bound 
(|47p for imaginary parts of resonances in Fig. [T^) 



large x and can be proved from semiclassical approxima- 
tion to ([57)1 . Physically, it is natural that all the states 
have a leakage smaller than a particle moving along the 
diameter of the circle. 

As was mentioned above, the resonances in Figs. [T^) 
and b) are not the only existing quasi-bound states of 
the dielectric disk. It is known (see e.g. [l^]) that func- 
tions Sfyi{x) have another series of zeros corresponding 
to the shape (or outer) resonances, im.p, also depend- 
ing on 2 quantum numbers m and p. At Fig. [2] some of 
these resonances are plotted. Notice the difference in the 
vertical scale of this figure in comparison with Figs. [T^) 
and b) . All resonances at Fig. [1^) are situated above the 
dashed line in Fig. [2l In [16] it was shown that positions 
of these outer resonances are described asymptotically by 
the equation 
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(48) 



Its solutions with the smallest modulus of the imaginary 
part have Xm,p ~ "m and are close to complex zeros of the 
Hankel function in the denominator of this equation. In 
|16| it is shown that in this case 

0{linZjn,p/ReZm,p) 



ri? — m? / z^ 



m,p 

(49) 
where Zm,p is a complex zero of Hm (z) function. It is 

known [30| that at the principal Riemann sheet oi Hm (z) 
cut along the negative part of the real axis such zeros have 
negative imaginary parts, lie symmetrically with respect 
to the imaginary axis, and for integer m there is [m/2] 
zeros with positive real part. These Hankel zeros with the 
lowest imaginary part can asymptotically be computed 
by the expansion (see e.g. [10 |) 



^m^p 



%lY 



1/3 



-27ri/3 



0{ill/m 



l/3^ 



(50) 



and r}p is the modulus of the p**^ zero of the Airy function 
(p = 1, 2, . . .). With a good precision these Airy zeros are 



described by the semiclassical formula 



and 



3/2 



7r(p 



b 



(51) 



which works well even for small p. In general, shape res- 
onances of dielectric circular cavity are situated at 0(1) 
distances from shape resonances of circular cavity with 
Dirichlet boundary conditions. 

From the above arguments it follows that the addi- 
tional resonances have quite large imaginary parts and 
are located below the curve 



Ima; < const IRexl"'^'^ 



(52) 



so they are well separated from the main resonances 
obeying (j46| . This bound has been proved [32], [33] for 
the scattering on any smooth (C°°) convex cavity pro- 
vided no trapping orbits exist. If the cavity boundary is 
not smooth only the weaker bound is known [341] 



— Ima; < const log ]Re: 



V. WEYL'S LAW 



(53) 



The circular cavity is useful to find the Weyl law ([28]) 
for the smooth part of the counting function. As usual 
(see e.g. j35i]), it can be extracted by considering pure 
imaginary values of momentum k — is and calculating 
the asymptotics of the Bessel functions from the known 
formulas (see e.g. [301] V In this manner one gets 



ISr. 



(is) 



2(m2-|-n2s2) 2(m2-)-s2) 



+ 1 
0{s-') 



(54) 



Changing the summation over m in (|40|) into the integra- 
tion from — oo to oo one obtains that two main terms of 
the smooth part of the right hand side of (dOj) have the 
following form 



d 



'dE 



L 



■Indet S(£;)> = ^{n^-l) + [r{n)-l]^+0{k-^) 

(55) 
where r{n) is defined in ([^5]) . 

But this formula gives the asymptotic behavior of the 
determinant of the full S-matrix which is not necessar- 
ily related directly with cavity resonances. According to 
P^ the S-matrix for the dielectric circular cavity is the 
product of 2 factors 



S^(i?) = SW(i^)SW(i?) 



where 



^^:i\E) 



Hil\x) 



(56) 



(57) 
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(58) 



It is easy to check that the S^"^ -matrix is the S-matrix for 
the scattering on a circular disk with the Dirichlet bound- 
ary conditions whose asymptotic behavior is known (see 
e.g. [33 and references therein) 
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The matrix element Sm (E) is the ratio of 2 functions: 
Sm{x) defined in (j37p and its complex conjugate. It is 
convenient to rewrite it as the ratio of two other functions 






where 



nJ!^{nx) — Jm{nx) 



Hl^^'jx) 
Hil\x) 



(60) 



(61) 



The factor x™^^ is introduced to insure that Fm{x) tends 
to a constant independent on x when m —^ 00. 

The function F.m{x) has 2 groups of zeros (cf. ([55]) ). 
The first includes all usual resonances which obeys the 
inequality (|46p and the second consists of additional res- 
onances (shape resonances) with large imaginary part 
([5^ . But it has also a series of poles coming from ze- 
ros of the Hankel function Hm (x). According to ([49]) 
the position of shape resonances in two leading orders 
differs from the corresponding zero of the Hankel func- 
tion only by a finite shift. From the bound ([5^ it follows 
also that for a finite 7 > there is only a finite number 
of additional resonances whose imaginary part obeys the 
inequality Imfc < — 7. 

These two arguments demonstrate that with a big pre- 
cision additional resonances and poles of the Hankel func- 
tions corresponding to shape resonances of the external 
Dirichlet problem cancel each other. It means that func- 
tion Fm{x) can be considered as having effectively only 
resonances with small imaginary parts. Consequently, 
the mean value of (Indet S'""^')' gives the average value 
of the counting function of resonances. From ([55]) it fol- 
lows that the later equals to the difference between ([55|) 
and ([591 



(^ Indet S(i)(i?)) 

(^ Indet S(i?))-(^ Indet &^'\E)) 



(62) 



which finally leads to the asymptotics ([?7|) . A careful sep- 
aration of the internal and external resonances in higher 
order corrections remains an open problem. 

The fact that in order to find the smooth part of the 
resonance counting function for a dielectric cavity in the 



leading order one has to divide the determinant of the 
S-matrix for the scattering on this cavity by the deter- 
minant of the S-matrix for the scattering on the same 
cavity but with the Dirichlct boundary condition may 
physically be argumcnted as follows. In the semiclassical 
limit a small part of a smooth boundary can be consid- 
ered as a straight segment. Resonances inside the cavity 
can be determined from the knowledge of the reflection 
coefficient from this part of the boundary as it is done 
in the Appendix [X] It is easy to see that the S-matrix 
corresponding to the reflection from outside the cavity 
differs just by a sign from the reflection coefficient inside 
the cavity. But the S-matrix equals —1 is exactly the S- 
matrix for the scattering on the cavity with the Dirichlet 
boundary conditions which explains the above statement. 
For a symmetric cavity it is convenient to split reso- 
nances according to their symmetry representations. Of- 
ten it equivalent to consider a smaller cavity where cer- 
tain boundaries are true dielectric boundaries but on oth- 
ers parts of the boundary one may have in the simplest 
cases either Dirichlet or Neumann boundary conditions. 
Now the total boundary contribution to the mean count- 
ing function can be written as follows 



-— (n(ijv - Ld) + r{n)Lo) k 
in 



(63) 



Here L^ and Ljj are the lengths of the boundary parts 
with respectively Neumann and Dirichlet boundary con- 
ditions and Lq is the length of the true dielectric inter- 
face. 

The knowledge of the spectrum for dielectric disk (cf. 
Fig. [1]) permits us to check the derived formula for the 
mean counting function. First we count how many res- 
onances exist with real part less than a fixed value, Re 
km < k. These resonances have to be taken into ac- 
count with their multiplicities. For the circular cavity 
states with m ^ are double degenerated and states 
with m = are simple. Then the resulting curves are fit 
by a polynomial of second degree with the highest term 
equals n'^Ak^/4TT. As for i? = 1, A = tt and L = 2tt, one 
gets the following fits: 



• For n — 1.5: 



N{k) = 



2.25. 



.5139fc + .1518 



• For n — 2: 

N{k) = k^ + .5446/c + .034 



(64) 



(65) 



The quality of the fits can be seen from Fig. [3] where 
the differences between numerically calculated counting 
functions and the above fits are plotted. The theoretical 
prediction for the coefficient of the linear in k term is 
f(n)/2 where f{n) is defined in ((28)) . Numerically one 
finds that these coefficients should be equal to .5125 for 
n = 1.5 and to .5420 for n = 2 in a good agreement with 
the above numerical fits. 
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FIG. 3: Difference between the exact counting functions 
of resonances of circular cavity with unit radius and their 
quadratic fits (|64p and ((65]): (a) for n = 1.5, (b) for n = 2. 



VI. PERIODIC ORBIT CONTRIBUTION 

To get oscillating contributions to the trace formula 
for a circular cavity we start from the expression (|40p . 
Expressing the Bessel function in ((37]) through the Han- 
kel functions of the first and the second kinds Jm{x) = 

{Hin (x) + Hm {x))/2 and performing straightforward 
transformations of s,„(a;) in ((57]) one gets a formal se- 
ries for the density of states 



d{E) - doiE) 
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E i^^^ 
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Bm{x) 



(66) 
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Using the Poisson summation formula, the sum over m 
can be substituted by the integration over m and the 
summation over a conjugate integer Af 



d{E) - do{E) 
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2mMm 



(67) 
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P„,{x)[R'^}{x)YEl,{x) 



S„i(.t) ttx ^ 
Here the following notations have been introduced 
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B„iix) = n — 75T-(nx) 



H, 



(1)' 



H, 



(2) 



H-, 



tt(^) ' 
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(68) 



(69) 



(70) 



(71) 



and 

Pm{x) = (n^ - 1) (^H^\nx)H^^\nx)Ara{x)B^{x)y^ . 

No approximation has been done in Eq. (j67p and it can 
serve as a starting point for derivations of dominant terms 
and corrections to them. In the semiclassical hmit x -^ 
oo it can be simpHfied by using standard asymptotics of 
the Hankel function (SOl 



H^^\^) 



(x2 - to2)1/4 



J<Prm{x) 
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bi{x) 



V: 



x^ — m^ 



+ 0(x-2) 



(73) 



where 

4>m{x) = \J rp-x'^ — ra^ 
and 



m arccos 



e-i (-' 



hi{x) 



1 
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(1 - xVm^) 



i\-\ 



(75) 



In this manner one concludes that, when x -^ co and 

\m\ < nx, E,n{x) -^ g2i0,„(a:)^ Pto(x) -^ \J Tp- — m? / x"^ , 

and 



Rt\x)-^R^{x) 



\/v?- — m? / x'^ — ^J\ — m? / x'^ 



sjn} — ni^ jx^ + -^/l — 777,2 /a 



(76) 

From these equations it foUows that the dominant peri- 
odic orbit contribution is 



B? 



CXD OO 



^'-'^^)-^. E E 



(77) 



A/=-oo r=\ 



dm\ n^ - _-i?;;^(a;)e'^"--('") + ^ 



where the action (the phase of the exponent) 

SM,r{m) = 2TrMm + 2r(j)„i{x) . (78) 

The last step, as usual, consists in the computation of the 
integral over ttz by the stationary phase method. The 
saddle point value of 777 = msp is determined from the 
equation 



_d_ 

dm 



SM.rim) = 



whose solution 



with 



iTT-sp = nx cos 6m r 



M 



IMa 



(79) 



(80) 



(81) 



corresponds geometrically to the periodic orbit of the cir- 
cular cavity in the shape of regular polygon with r ver- 
tices going around the center M < r times. If M and 
r are co-prime integers the periodic orbit is primitive. 
Otherwise, it corresponds to a primitive periodic orbit 
repeated d times where d is the largest common factor of 
M and r. 

The action ((78)) can be expanded into a series of devi- 
ation from from the saddle point value 777 — rrisp + 5m. 
One gets 

SmA^sp + 5m) = nklp - ^r + ^Sm"^ + 0{5m^) (82) 

where Ip is the periodic orbit length 



Ip = 2rRsiii9M,r 



(83) 



and (72 is the second derivative of the action computed at 
the saddle point 



2r 



92 



nx sm Um r 



(84) 



Computing the integral over 777 by the saddle point ap- 
proximation one gets the expected result (|25p 



d(o-)(£;) = !^y^_^_i?;e'["'='''-''"/2+-/4] ^.c.c. 



^ tE^ y/^T^klp 



where 



Ap — ttR^ sin^ 6m., 



(85) 



(86) 



is the area occupied by the considered periodic orbit fam- 
ily and Rp is the Fresnel reflection coefficient ([S|) calcu- 
lated at the incidence angle 



77 — '^M,r 



(87) 



As classical motion inside a circle conserves the angle 
with the boundary, the mean value in (|25p equals the 
product of reflection coefficients in all incident points. 

At Fig. 3^) the length density (i.e. the Fourier trans- 
form of the trace formula pTjl ) is presented for the circu- 
lar dielectric cavity with the refraction coefficient equal 
to 1.5. The vertical solid lines indicate the lengths of 
simplest periodic orbits. From the left to the right these 
lines correspond to, respectively, triangle, square, pen- 
tagon, hexagon, heptagon, and octagon. For clarity the 
dotted line at this flgure represents the length density 
for the closed circular cavity with the Dirichlet bound- 
ary conditions. As expected the both models have peaks 
at periodic orbit lengths but their amplitudes are differ- 
ent. In particular, as for 77 — 1.5 the critical angle is 
of the order of 42°, triangular periodic orbit is not con- 
flned and the corresponding peak for dielectric cavity is 
practically invisible. Others peaks have similar heights. 
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a) b) 

FIG. 4: (a) Solid line: length density for dielectric circular 
cavity with n — 1.5 and R = 1. Dotted line: the same but 
for circular cavity with the Dirichlet boundary conditions, (b) 
Goos-Hanchen shift associated with the square periodic orbit. 
Solid line corresponds to the dielectric disk and the dotted 
line represents the circular closed cavity. Small vertical line 
indicates the maximum of length spectrum for dielectric disk. 



CORRECTIONS TO ASYMPTOTIC 
REGIME 



VII. 



According to p5l) the contribution of a periodic orbit 
for an open dielectric cavity differs from a closed cavity 
only by the value of the total reflection coefficient. Thus, 
when all incidence angles for a periodic orbit are greater 
than the critical angle ([9]), the modulus of its contribution 
to the trace formula has to be the same as for the closed 
cavity. The curves at Fig.[¥]have been normalized in such 
a way that their peaks at pentagonal periodic orbit are 
the same heights. But then the amplitudes corresponding 
to the square orbit are clearly different (dielectric cavity 
peak is of the order of 0.7 of the close cavity peak) which 
contradicts the above prediction. 

Other important corrections are the so-called Goos- 
Hanchen shift [lj|, [3a, [331 which manifests itself as a 
slight difference of a visible length from its geometric 
length for an orbit reflected from a dielectric interface 
and the Fresnel filtering [38| , which manifests itself as a 
slight shift of the angle of emission near the critical angle 
for a dieletric boundary. At Fig. UJj) a region of Fig. 3^) 
close to the maximum of the length density around the 
square periodic orbit is magnified. The peak of the close 
cavity represented by dotted line is, as expected, centered 
at the geometric length of this orbit (/ — 4-\/2i? ~ 5.6568 
for R = 1). But the peak associated with the dielectric 
circular cavity is clearly shifted to the right. Numeri- 
cally, this shift is small and for the considered interval of 
momentum it is of the order of 0.004i?. 

The reason of these discrepancies is related with the 
fact that asymptotics of the Hankel function ([75)1 is not 
valid when m is close to x (see e.g. ,3Q|]). The uniform 
approximation for large m and all x is more complicated 
and is given e.g. by the Langer formula [30| 



HilH^)-[l 



arctan to 



1/2 



H[%{z)e'-/' 



0{m 



-4/3\ 



(88) 



where z = mioj — arctan uj) and uj = ^ x"^ jm? — 1. 

Non- uniform character of this expression is clearly seen 
from the following asymptotics useful for further estima- 
tions 
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(89) 
(90) 
(91) 



.4592-^.79541. 



where a = 2-2/335/6^-ip2(2/3)e2-i/3 

To calculate carefully higher order corrections to the 
trace formula ([55)1 is not an easy task. The main diffi- 
culty lies in the fact that the stationary phase method 
on which trace formulas are based on cannot, strictly 
speaking, be applied when a periodic orbit has incident 
angles close to the critical one. The point is that in the 
stationary phase method it is implicitly assumed that 
non-exponential terms are changed much slowly than the 
action in the exponent. This is the case for the usual bil- 
liard systems where the whole perturbation series for the 
trace formula has been constructed [33] , [40] ■ 

But as follows from ([55)1 for dielectric cavities exact 
refiection coefficient ([M)) near the critical angle changes 
at a scale of 6rh ^ x^^^ but the action changes at a 



different scale of the order of 5ifi 



r-l/2 



When X — > 00 the former is much smaller than the 
latter and the stationary phase method when a trajec- 
tory hits a boundary with angle close to the critical one 
can not be justified. To avoid this difficulty one can use 
the Langer approximation ([55]) in the transitional region 
and tabulate the integrals without other approximations. 
We will not perform these computations here. Instead, 
we compare the saddle point result with numerically cal- 
culated integral for the square and pentagonal periodic 
orbits. For refractive index n = 1.5 incident angles of 
the square periodic orbit are close to the critical angle 
but for the pentagonal orbit they are far apart and the 
difference of these two cases will give the accuracy of the 
pure saddle point results. 

From ((67|) it follows that the full contribution of a pe- 
riodic orbit with fixed numbers M and r is given by the 
integral 

(92) 



dp{E) = S'^^^^-'^/^+^'^^JM^rix) 



where 



JmA^) 



R^ 



P^{x)[Rt\x)rEl^{x) 



(93) 



^g-i[nfcip-r7r/2+7r/4]^^ 



In the pure stationary phase approximation integral (|93l) 
equals 

j):F}i^) - -^Pra,M[Rtl{x)]\r^ (94) 

' TT^a;^ ^ V <?2 
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which in the semiclassical hmit tends to 



^M^r \^) 



^T^y2^M 



^BL 



V 



(95) 



p 



as in ((85|) . 

At the first sight it seems that a better approxima- 
tion consists in the calculation of the prefactor ((72|l and 
the reflection coefficient (|69| in ([94]) directly at the saddle 
point ([80| without the asymptotic formula ([73|) . Roughly 
speaking, this is equivalent to use for the reflection on 
a curved surface instead of the Fresnel diffraction coeffi- 
cient ([U valid for an infinite plane interface the reflection 
coefficient (|69|) computed at the saddle point value of the 
incidence angle 9 as in ([87| . When this angle is close to 
the critical angle from ([M|) one gets 



Rrn{0,x) 



(1)'/ 



rWt 



ncos9 + iHii>'{x)/H^'{x) 
n cos 61 - iHil^'{x)/H^\x) 



(96) 



where the value of the angular momenta m — m{9) is 
related with angle 6 as it follows from the saddle point 
relation §0^ 



m{6) — nxsiiiO 



(97) 



This type of modification first has been proposed in 37| . 

Nevertheless, we found that this approximation numer- 
ically is not better than the simple semiclassical formula 
([95)1 . This is because the effective reflection coefficient 
([96]) changes too quickly close to critical angle. There- 
fore, when a region which gives a dominant contribution 
to the integral ((93)) is close but not equal to the saddle 
point in the exponent, expression ([M)) differs consider- 
ably from the correct one. For this reason we prefer to 
compare numerical values of the integral (|93p not with 
([94)1 but with simpler expression ((95)) . 

At Fig. [5] the results of numerical calculations of the 
integral (|93|1 for square and pentagonal periodic orbits 
are presented. At Fig. [5^) the square of the modulus of 
this integral divided by the semiclassical expression (|95)) 
is given. As expected, contribution of the pentagonal 
periodic orbit is well described by semiclassical formula. 
But the contribution of the square periodic orbit differs 
considerably from the asymptotic result which explains 
quantitatively the observed differences of the square orbit 
heights for closed and dielectric cavities seen at Fig. |4| 

The computed corrections permit also to explain the 
analog of the Goos-Hanchen shift observed at Fig. [JJd) . 
Let us write the periodic orbit contribution in the form 



ipiE) = \dp{E)[ 



,i[fcnip-rTr/2+7r/4+*p(x)] 



(98) 



where the phase ^p(a;) comes from the integral ()93p . 

To compute the length spectrum one has to calculate 
the Fourier harmonics of this expression (|32)) . The shift 



of the peak position, Al = I — Ip can be estimated from 
the saddle point in the exponent 



Al _ d-^pix) 
^"R ~ dx 



(99) 
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FIG. 5: (a) Modulus square of the integral (I93|l divided by its 
semiclassical expression (|95|l . Solid line indicates the contri- 
bution of the square periodic orbit. Dashed line is the contri- 
bution of pentagonal orbit, (b) The derivatives of the phase 
of the same integral: solid and dashed lines correspond re- 
spectively to the square and pentagonal orbits. 



At Fig. [5)d) the value of this derivative computed numeri- 
cally from the knowledge of ([M)) is plotted for the square 
and pentagonal orbits. For the square orbit its value is 
larger than for pentagonal one and is of the same order 
as the observed shift in Fig. Hja). As we are looking for 
a very small effect, mutual influence of different periodic 
orbits has to be taken into account when a precise deter- 
mination of the Goos-Hanchen shift is required. 



VIII. CONCLUSION 

The goal of the paper is the construction of the trace 
formula for open dielectric cavities. By using the Krein 
formula it was demonstrated that such trace formula can 
be written in the form ([50)1 



1 



n ^ {k-ReE„ 



Im E„ 



{ImE„ 



d{E)+Y,{dp{E)+d;{E)) 



(100) 



where in the left-hand side the sum is taken over all inter- 
nal resonances and in the right-hand side the summation 
includes all periodic orbits for the free motion inside the 
cavity. The term d{E) is the mean density of these res- 
onances given by the derivative of ()27p . Its two main 
terms are 



diE) 



An' 
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(101) 



where k = VE, A and L are the area and the perimeter 
of the cavity, and f{n) is as follows 



n di „, , 1 / 



dt 



i2 + l 
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where R{t) is the reflection coefficient at imaginary mo- 
mentum defined in (|29p . Terms related with R{t) cor- 
respond to the reflections from inside and outside the 
boundary but 1 in (|102p appears due to the necessity of 
removing the S-matrix for the scattering on the cavity 
with the Dirichlet boundary conditions. 

The contribution from a particular periodic orbit, 
dp{E), to the trace formula (|100p differs from the same 
closed billiard only by the product over Fresnel reflection 
coefficients ^ computed at all reflection points as in ipi]) 
for isolated orbits or their average over a periodic orbit 
family as in (P5|) . 

We have discussed also corrections to the above asymp- 
totic semiclassical formulas. Due to rapid changes of the 
reflection coefficient close to the critical angle corrections 
for dielectric cavities are larger than for closed cavities. 
In particular, higher order corrections permit to estimate 
the analog of the Goos-Hanchen shift for the peak posi- 
tion in dielectric cavity length density. 

Obtained formulas were compared with numerical cal- 
culation of resonances for the circular dielectric cavity 
and in all cases a good agreement was observed. 
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APPENDIX A: KREIN FORMULA FOR ONE 
DIMENSIONAL SYSTEMS 

To illustrate the Krein formula in the simplest case 
of a one-dimensional system consider a particle moving 
along the positive half-axis in a potential without bound 
states. For simplicity, assume that (i) ip{Q) = and (ii) 
the potential is zero for distances larger than a. The cor- 
responding quantum problem has only continuous spec- 
trum. To get artificially the discrete spectrum let us im- 
pose that at some large distance R^ a the wave function 
obeys 



ip{R) = . 



(Al) 



Eigenfunctions of a such system at large distances x > a 
have the form of plane waves 



'ip{x) ex sm{kx + S{E)), 



(A2) 



where k = VE and the phase shift S{E) contains all 
information about the potential. 

Eigcnmomenta, k„i, of this system are determined from 
the condition (jAip and they read 



k^R + 5{E) 



nm 



(A3) 



FIG. 6: 1-dimensional dielectric cavity. Refractive indices are 
indicated above the axis. 



where integer m counts the number of states. By defini- 
tion, the density of states, d{E), is the number of states 
within an interval dE 



^ ' dE 2ttVE 



1 dSjE) 
^ dE 



(A4) 



But the function (JA2| with S{E) = is the eigenfunction 
of the free problem without the potential. Therefore, the 
first term in the last equation describes the free density 
of states 



doiE) 
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2TrVE 
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The excess density due to the presence of the potential is 
the difference between the level density with the potential 
and without it. From (IA4I) it follows that 



d{E) - do{E) = 



1 dSjE) 
n dE 



(A6) 



Instead of the phase shift it is instructive to introduce 
the S-matrix for the scattering on the potential defined 
as the ratio of the amplitude of the outgoing wave to the 
amplitude of the incoming wave. For the function (|A2p 
one has 



SiE) 



^2i5{E) 



Therefore, the excess density of states reads 



d{E) - do{E) 



1 d\nS{E) 
2^1 dE 



(A7) 



(A8) 



which is the 1-dimcnsional version of the Krein formula 

It is also straightforward to compute explicitly all nec- 
essary quantities for a 1-dimensional dielectric cavity [lO| . 
Consider a dielectric segment of size L with the refrac- 
tive index n in a media with the refractive index 1 as 
at Fig. [6l After imposing condition (jAip one gets that 
the density of states d{E) has to be determined from the 
formula 

1 r^ 

d{E)^ / lmn'^{x)G{x,x)dx (A9) 

■^ Jq 

where n{x) — n inside the cavity and n{x) = 1 outside 
it, and G(x, y) is the Green function of our problem [lO|. 
Using 



doiE)^-- f Im Gf\x,x)dx = -^ 



(AlO) 
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where G'''^(a;, y) = e'''l^~^'/2ifc is the free Green function 
after some calculations one finds that for antisymmetric 
states when i? — > cx) 



diE)-do{E) = - — 



nL 



1-H 



L 



27rfc 27rfc ll + re' 



2inkL\2 



(All) 



where r is the Fresnel reflection coefficient from dielectric 
boundary 



n-1 
n + 1 



(A12) 



Correspondingly, the S-matrix for antisymmetric states 
is 



SiE) 



^2i{n-l)kL 



1 



re 



-2inkL 



1 + re 



2mkL ' 



(A13) 



which, of course, agrees with 1-dimensional Krein for- 
mula (|A8|) . This S-matrix has an infinite number of 
poles (or resonances) corresponding to the solutions of 
the equation 1 + j-e^'"*^™^ = or explicitly 



riL^ 2' 2nL 



In- 



1 



(A14) 



where m is an integer. 
Using the identity 

l-r2 



1 + r^ + 2rcos(a;) 



E 

j=-oo 



i-r) 



liLya 



(A15) 



valid for |r| < 1 one can rewrite (jAlip as the trace for- 
mula for resonances in the considered simplest case 



^ oo 



ImA;„ 



TT •' — ' (k — Kekm)^ + (Imfc„ 

rn—oo ^ ^ ^ 



^1^ ,!:<-) 



j\^2ijnLk 



(A16) 



The term with j = in this equation corresponds to the 
mean density of resonances and terms with j y^ describe 
the contributions from the unique primitive periodic or- 
bit and its repetitions. 



Notice that the transition from the Krein formula 
(|A11[) to the trace formula (jA16p is done by removing 
from the S-matrix (PT3|l the factor So{E) = -e'^'''^ 
which is exactly the S-matrix for the scattering on the 
same cavity but with the Dirichlet boundary condition 
at the right end of the cavity x = L. This is the same 
phenomenon as for the circular dielectric cavity discussed 
in Section [3] 
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